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FOREWORD 

This  report  describes  theoretical  and  experimental  analyses  developed  by 
staff  of  the  University  of  Colorado  at  Boulder  in  collaboration  with  the 
Electromagnetic  Fields  Division  of  the  National  Bureau  of  Standards  (NBS), 
under  a  contract  sponsored  by  NBS.  Professor  David  C.  Chang  heads  the 
University  team.  Dr.  Mark  T.  Ma  of  NBS  serves  as  the  technical  contract 
monitor.  The  period  covered  by  this  report  extends  from  July  1980  to  July 
1981. 

The  work  presented  in  this  report  represents  a  further  aspect  of 
establishing  a  theoretical  basis  for  the  technical  analyses  of  transverse 
electromagnetic  (TEM)  transmission  line  cells.  The  general  purpose  of 
pursuing  theoretical  studies  is  to  evaluate  the  use  of  TEM  cells  for  (1) 
measuring  the  total  rf  radiated  power  by  a  device  inserted  into  the  cell  for 
test,  or  (2)  performing  necessary  susceptibility  tests  on  a  small  electronic 
device. 

The  particular  topic  addressed  herein  is  to  determine  the  input  impedance 
of  a  probe  antenna  which  may  be  inserted  into  a  TEM  cell  to  measure  or  excite 
fields.  The  formulation  of  the  problem  is  based  on  a  variational  principle, 
which  insures  that  the  input  impedance  so  obtained  is  stationary  for  small 
arbitrary  variations  in  the  probe  current  distribution  about  its  true  value. 
The  evaluation  of  the  variational  integral  requires  a  knowledge  of  the  field 
distribution  inside  the  cell  excited  by  a  vertical  short  electric  dipole. 
This  latter  problem  was  the  subject  of  a  previous  report  by  the  same  authors 
[4],  the  results  of  which  are  readily  adapted  to  the  present  problem.  Thus, 
the  input  impedance  treated  in  this  report  represents  an  extension  to  the 
previous  studies. 

The  analysis  of  the  present  problem  is  significantly  simplified  by 
assuming  that  the  gap  between  the  inner  septum  and  the  outer  walls  of  the  cell 
is  electrically  small,  and  that  the  operating  frequency  is  such  that  no  higher 
order  modes  are  excited  inside  the  cell.  These  assumptions  are  consistent 
with  the  practical  cell  utilization  and  design  principle. 

The  final  impedance  expression  is  shown  to  consist  of  two  distinct  terms, 
which  can  be  clearly  identified  as  the  respective  contributions  by  the 
ordinary  rectangular  waveguide  and  the  gap  perturbation. 
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INPUT  IMPEDANCE  OF  A  PROBE  ANTENNA  EXCITING  A  TEM  CELL 

P.  F.  Wilson,  D.  C.  Chang  and  M.  T.  Ma 

The  input  impedance  of  a  probe  antenna  exciting  a  transverse 
electromagnetic  (TEM)  transmission  line  cell  is  formulated  by  a 
variational  approach.  The  formulation  also  utilizes  the  results 
from  a  previous  work  on  the  field  distribution  inside  a  TEM  cell 
excited  by  a  vertical  electrical  Hertzian  dipole.  The  final 
result  of  imoedance  is  shown  to  consist  of  two  distinct  terms, 
which  are  respectively  contributed  by  the  ordinary  rectangular 
waveguide  and  the  gap  perturbation.  Numerical  results  for  both 
the  real  and  imaginary  parts  of  the  impedance  are  given.  The 
resistive  part  is  found  to  be  proportional  to  the  square  of  the 
probe  length,  and  the  reactive  part  largely  capacitive. 

Key  words:  Green's  function;  input  impedance;  probe  antenna; 
radiation  resistance;  rectangular  coaxial  transmission  line;  TEM 
cell;  variational  method. 

1.   INTRODUCTION 


The  National  Bureau  of  Standards  (NBS)  has  been  interested  in  developing  a 
method  for  testing  the  emission  and  susceptibility  properties  of  electronic 
equipment.  One  approach  to  this  problem  is  the  transverse  electromagnetic 
(TEM)  transmission  line  cell.  The  TEM  cell,  shown  in  figure  1,  consists  of  a 
section  of  rectangular  coaxial  tramission  line  (RCTL)  coupled  at  each  end  to 
standard  50ft  cylindrical  coaxial  line  via  a  tapered  section.  The  RCTL  section 
supports  a  dominant  TEM  mode  and  provides  an  isolated,  standard  test  field 
environment.  The  TEM  cell  results  may  be  related  to  the  free  space 
environment  via  results  given  by  Tippet  [1]  and  Sreenivasiah  [2], 

This  report  concerns  the  manner  in  which  an  idealized  coaxial  probe 
excites  an  infinite  section  of  RCTL.  The  results  may  be  used  to  study  the 
effect  of  coaxial  probes  inserted  into  the  TEM  cell  either  to  excite,  or  to 
measure  fields.  In  addition,  the  configuration  studied  may  be  useful  in 
modeling  feed  lines  to  electronic  equipment.  The  method  of  solution  is  to 
formulate  a  variational  expression  for  the  input  impedance  of  the  probe  at  the 
aperture  plane,  as  was  done  by  Collin  [3]  in  analyzing  the  coaxial  probe 
excitation  of  a  rectangular  waveguide.  The  evaluation  of  the  variational 
integral  requires  that  the  Green's  function  corresponding  to  the  fields 
excited  in  an  RCTL  by  a  vertical  electric  dipole  (VED)  be  known.  The  VED 
excitation  of  an  RCTL  was  the  subject  of  a  previous  report  by  the  authors  [4], 
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the  results  of  which  will  be  adapted  to  the  present  problem.  The  analysis  of 
the  Green's  function  is  significantly  simplified  by  assuming  that  the  gap 
between  the  inner  septum  and  the  outer  walls,  is  electrically  small,  and  that 
the  operating  frequency  is  not  too  large.  Both  these  premises  are  consistent 
with  the  TEM  cell  usage. 

Evaluating  the  variational  impedance  integral  results  in  two  physically 
distinct  terms,  which  may  be  referred  to  as  the  ordinary  rectangular  waveguide 
contribution,  and  the  gap  perturbation  contribution.  The  ordinary  term  gives 
the  input  reactance  due  to  a  coaxial  probe  in  a  rectangular  waveguide  with  the 
same  dimensions  as  that  of  the  RCTL's  upper  chamber.  Since  the  dimensions  of 
the  cell  are  such  that  the  rectangular  waveguide  modes  in  the  upper  chamber 
would  be  cutoff,  the  result  is  a  highly  capacitive  term  for  which  the  analysis 
yields  an  approximate  analytic  expression.  The  gap  perturbation  term 
represents  the  change  from  the  ordinary  result  due  to  the  presence  of  the  gap 
fields,  and  contains  the  radiation  resistance  resulting  from  the  dominant  TEM 
mode.  The  perturbation  term,  as  derived,  remains  in  an  integral  form  which 
must  be  treated  numerically.  However,  if  the  TEM  mode  radiation  resistance  is 
extracted,  the  remaining  reactive  integral  is  negligible  when  compared  to  the 
highly  capacitive  ordinary  contribution.  Thus  neglecting  the  perturbation 
reactance  results  in  a  computationally  simple  algebraic  form  for  the  probe 
input  impedance. 

Section  2  of  this  report  will  define  the  variational  impedance  integral  to 
be  evaluated  and  discuss  the  necessary  probe  idealizations.  The  Green's 
function  associated  with  the  VED  excitation  of  an  RCTL  is  examined  in  Section 
3.  Section  4  gives  the  result  of  the  evaluation  of  the  impedance  integral  and 
examines  the  distinct  physical  nature  of  the  resulting  terms.  Certain  special 
cases  are  used  to  verify  the  solution.  These  checks  are  examined  in  Section  5 
along  with  numerical  results.  Many  details  of  the  analysis  and  a  listing  of 
the  FORTRAN  program  written  to  evaluate  the  impedance  terms  are  left  to  the 
Appendices. 


2.  VARIATIONAL  FORMULA  FOR  THE  PROBE  INPUT  IMPEDANCE 

The  RCTL  cross  section  and  the  coaxial-line  probe  antenna  are  depicted  in 
figure  2.  The  septum,  of  width  2w,  is  allowed  to  be  vertically  offset  with 
the  heights  of  the  upper  and  lower  chambers  given  by  b-.,  and  bp, 
respectively.  The  width  of  the  upper  chamber  is  2a-i,  and  the  width  of  the 
lower  chamber  is  2ap.  The  gap  between  the  septum  and  the  outer  wall  is  g. 
The  probe,  of  length  d  and  radius  t,  is  located  in  the  upper  chamber.  A  TEM 
cell  would  have  an  RCTL  cross  section  with  equal  chamber  widths  and  previous 
discussions  have  assumed  this  geometry.  However,  in  order  to  avoid  any 
degeneracies  caused  by  letting  a-i  =  ao  =  a,  the  analysis  will  initially  assume 
general  chamber  widths.  The  time  convention  will  be  exp(-icot). 

If  the  coaxial -line  aperture  is  small,  it  is  reasonable  to  neglect  the 
higher  order  modes  excited  at  the  aperture  discontinuity  and  assume  only 
incident  and  reflected  TEM  modes  exist  in  the  coaxial-line.  As  shown  by 
Collin  [3]  in  his  analysis  of  the  coaxial-line  probe  antenna  excitation  of  an 
ordinary  rectangular  waveguide,  the  above  approximation  allows  the  antenna 
input  impedance  Z  to  be  expressed  as  follows 

z  =  4"  /  /  JfxJ-KCx.x'-O'JCx'Jds'ds         (1) 
I   S  S' 

where  x  and  x'  denote  source  and  observation  points,  J(x)  is  the  unknown 
current  density  on  the  probe,  K(x,x')  is  a  dyadic  Green's  function 
corresponding  to  the  electric  field  excited  by  a  unit  strength  dipole  source, 
and  It  is  the  total  current  at  the  aperture  plane.  S  denotes  the  probe 
surface  plus  the  aperture  plane  which  carries  an  equivalent  electric  current 

J(x)  =  nxH(x)  where  n  is  a  unit  vector  in  the  outward  normal  direction. 

Experimental  studies  have  shown  that  Z  is  effectively  independent  of  the 
aperture  size  for  electrically  small  apertures  [3].  Therefore  it  is  reasonable 
to  idealize  the  probe  as  driven  by  a  y-directed  E-field  concentrated  at  the 
base.  Thus  we  need  consider  a  y-directed  current  only,  which  reduces 
J(x)  to  a  J(x),  and  K(x,x')  to  ayK(x,x')ay  with  ay  being  the  y-directed  unit 
vector. 


The  above  impedance  expression  may  be  shown  to  be  variational  in  J(x). 
That  is,  a  first  order  approximation  for  J(x)  leads  to  a  second  order 
approximation  for  Z.  At  normal  TEM  cell  operating  frequencies,  the  length  of 
the  probe  will  be  small  compared  to  a  wavelength.  Therefore,  it  is  reasonable 
to  assume  that  the  current  distribution  is  approximately  sinusoidal, 

I 
J(x)  = -^t  sin  k(brd-y)6(P-t)  (2) 

where  the  delta  function  is  interpreted  in  the  sense  of  a  generalized 
function.  The  form  of  the  Green's  function  K(x,x')  will  be  considered  next. 

3.  THE  GREEN'S  FUNCTION  K(x,x') 

The  appropriate  Green's  function  for  the  special  case  of  equal  chamber 
widths  has  been  derived  previously  [4].  However,  the  modifications  due  to 
generalizing  the  chamber  widths  are  not  difficult.  Since  a  VED  in  an  RCTL 
will  excite  both  transverse  electric  (TE  or  h-type),  and  transverse  magnetic 
(TM  or  e-type)  modes,  both  cases  need  to  be  treated.  For  each  modal  type  we 
begin  by  formulating  an  integral  equation  for  the  unknown  field  in  the  gap. 
These  integral  equations  could  be  solved  by  Chebyshev  polynomial  expansion 
methods  [5],  However,  a  less  cumbersome  solution  may  be  achieved  by  making 
the  following  pair  of  assumptions.  First,  the  gap  is  assumed  to  be  small,  and 
second,  the  frequency  is  not  allowed  to  be  too  high.  These  assumptions  lead 
to  an  approximate  integral  equation  which  may  be  solved  directly.  The  small 
gap  assumption  amounts  to  discarding  terms  which  are  of  the  order  G  JinG,  where 
G  is  a  normalized  gap  parameter  defined  by  G=7Tg/2a,  a=min(a-^,a2) .  Up  to  the 
order  of  this  approximation  only  the  TE  type  modes  contribute  significantly  to 
the  gap  field. 

The  derivation  of  the  Green's  function  begins  with  the  removal  of  the  z- 
dependence  via  a  Fourier-transform  defined  by 


K(j)(x,x')  =  —  /  K(j)(x"  ,x')elaZda         (3) 

2tt  t   t 


where  j=l ,2  refer  to  the  upper,  and  lower  chambers  respectively,   the  subscript  t 
denotes  the  transverse  plane,  a  is  the  transform  variable  representing  a 
continuum  of  propagation  factors,   and  the  superscript  tilde  signifies  a  function 
in  the  transform  domain.     The  analysis  outlined  above  yields  the  following 
result 


where  K       represents  the  Green's  function  expected  if  there  were  no  gap,   i.e., 
the  ordinary  rectangular  waveguide  result,   and  i<2     gives  the  gap  field 
perturbation  contribution.     K         and  K2       derived  in  Appendix  A  are  found  to  be 

2     2 
2  A   [M   +a  ] 

~(j)    -      -  ia»M      a  -      -  n      j  (j)    -      - 

K  x    ,x'      =   -6       —  6  x   -x'      +       I       g         x    ,x') 

1         t     t  lj       2    l   2         t     t  L         2         2     2     mn       t     t  J 

k       c  m,n=0  M     +  N   -?  (5) 

j         j 


^»(x    x^-J^U.Jf  (x)f  (x;>  (6) 

2aiY 

where  u>  =  2irf  is  the  angular  frequency,  p  is  the  permeability,  k  =  ui(\ie)'2  is 
the  TEM  mode  propagation  constant,  6^  is  the  Kronecker  delta  function,  rj-  - 
k2-a  ,  M-  =  nrn/2aj,   and  Nj  =  nir/bj.     The  remaining  functions  are  defined  by 


(i)   -     -  2 

9       (xtSx!)   =   (— — )sinM.(x+a.)sinM.(x'+a.)cosN.ycosN.y'  (7) 

mntt  a.b.  jj  J  J  J  J 

(      8/a, a,              2       ,             cotK       b.         ,    ) 
L(o,   .L—Ll    tI    I    1    I      %+  +  ■+[  (8) 

o  m 


M.cosK  J    (b.-y) 

f.(xj   =  I       J,.,  m     , ."? sin^-cosM.xJ     (M.g)  (9) 

m       K       sinK       b. 
on  m       j 

Kij)     =  ?2  "  Mj  (10) 

with  m0  denoting  the  summation  over  odd  m  only   (m=l ,3,5. . .) ,   and  J     is  the 
zeroth  order  Bessel    function.     Again  the  delta  function  should  be  understood 
in  the  sense  of  a  generalized  function. 

4.      THE    INPUT   IMPEDANCE   Z 

4.1.     The  Green's  function  terms  K,J     and  k!j     will    lead  to  two  distinct 
impedance  contributions  Zj  and  Z2  which  represent  the  ordinary  and 
perturbation  contributions  to  the  probe  input  impedance.     The  total    impedance 
is  the  sum  of  these  two  terms.     Combining  the  form  of  the  Green's  function 
given  by   (4),  with  the  impedance  integral    (1),  yields  Z  =  Z^  +  1^,  where 

00 
Z     =  — ^    /  /     J^)   /     k!1^   ,x')eloZ  J(x')    da  ds'    ds  (12) 

2ttI.      S   S' 


Z2  =-^-    /  /    J(x)     /     K^^x^x^e^vKx')   da  ds'   ds  (13) 

2TrIt     S  S' 

The  integrals  defined  for  Z^  and  Z^  are  evaluated  in  Appendices  B  and  C  for  a 
centrally  located  probe.     The  approximate  form  found  for  Z^  depends  upon  two 
assumptions:     (1)   the  probe  is  taken  to  be  thin,   i.e.,     irt/2a.<<l,   and  (2)   the 
frequency  is  restricted  so  that  (ka^)2  is  small.     More  precisely,   terms  in  t 
of  the  order  (irt/2a1)2  are  discarded  and  terms  in  ka^  of  the  order  (k/2M1)4 


are  discarded.     The  evaluation  of  Z2  involves  no  approximations.     However,   Z7 
remains  in  an  integral    form.     The  results  of  Appendices  B  and  C  are  as  follows 

2  2 

Z     „  J3L_  tan  2J^£n   (-4-1   +^V  (4.207175) 
1       ^TTb^2  2!  rt  / 

9     »     ,        sin2N.d/2  \2  K   [t(N2-K2  ^    ) 
-2k2     I     (l-         ^  °     2   1     2  (14) 

n=l   \       sin     kd/2  /  fq  -  r         J 

2  »  1 

Z      =  ^^-  CSC2kd   /  -^^h    ,(o,d)    I      h    (a,d)J    [tfa^V2]    da                 (15) 

2,2  2           m                      mo               1 

4a,                   -«  c       m                    m 

1  00 
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M.    sin^f-J    (M.g)      (cos  kd  -  cos  K(1)d) 
i/i\i  ^01  m 

Va'd)   "  ,(1)      .     „(1)   . 7T—-2 

Km       S1  n  Km       bl  Ml  +  a 

where  KQ(x)   is  a  modified  Bessel    function. 

The  perturbation  term  1^  contains  the  radiation  resistance  due  to  the 
propagating  TEN!  mode,  which  may  be  extracted  by  evaluating  the  half  residues 
at     a  =  ±k  due  to  the  forward  and  backward  propagating  TEM  modes.     The  TEM 
mode  radiation  resistance  is   found  to  be: 

,2 

771-1  k  7  ?    ?  1/ 

R  =  -——  esc  kd     L(k)   I  h   ,(k,d)  I  h   (k,d)J   [t(k  +M  V2]  (17) 

2  m '      '       L     m     '       0 

4a,  m  m 

1  00 

where  ti     =   12(k. 


If  we  now  let  a^^^  s0  that  we  are  considering  a  normal   TEM  cell 

r 
by  [4] 


geometry,  and  recall  that  the  characteristic  impedance  2Q   of  an  RCTL  is  given 


n  it 
Zc  =-g-L(k)  (18) 


The  radiation  resistance  may  be  written 


?k  2  2     2  l/o 

R=^-Z     esc  kd  ^     h   ,(k,d)   \  h   (k,d)J   [t(k  +M  )' 2  ]  (19) 

2       c  m  mo 

a  mm 

o  o 

Because  only  the  TEM  mode  is  allowed  to  propagate,   R  represents  the  total 
resistance.     We  may  also  express  Z  as  the  sum  of  the  radiation   resistance  R 
and  a  reactance     i(X+6X),  with  X  and  6X  defined  by 


n            ~  .  i  e       <<  2,2 

o          2  kd  )     ,4a.  a  k     , .  „„_,,_,,_. 

X  =  ^-r-j-  tan    —  Un(— )  +  — — -  (4.207175) 

2ttd.  k               2  (        ut  2 


1 


TT 


9     -     /        sinVd/2  >2  K  [t(N2-k2)/2]  j 

-*  *  t -- r1-) -w- }  (2o) 

n=l    v        sin  kd/2     J         N     -  k  J 

2  «>  V2 

6X  =  ^_  csc2kd  P  /    t^-  J.     h   ,(o,d)  I     h   (a,d)J   [t(a2+M2)   ]     da 
.2  J         2  m  m  0 

4a  -°°       c      m  m 

00  (21) 


where  the  integral  is  now  understood  in  the  principal  value  sense.  The 
dominant  reactive  contribution  comes  from  the  sum  over  n  and  therefore  the 
ractance  is  highly  capacitive.  The  capacitive  behavior  is  expected  from 
comparisons  to  either  a  monopole  in  free  space,  or  a  similar  probe  in  a 
rectangular  guide.  For  typical  TEM  cell  geometries,  the  principal  value 
integral  does  not  significantly  add  to  the  reactance.  Therefore,  a  yery   good 
algebraic  approximation  for  the  reactance  is  simply  X  which  allows  us  to  write 
Z  as  follows: 

Z  *  R  +  iX,  (22) 


where  R  and  X  are  given  respectively  in  (19)  and  (20). 


5.  NUMERICAL  RESULTS  AND  DISCUSSIONS 

■ 

In  order  to  check  our  result  we  would  like  to  compare  our  expression  to 
certain  special  cases  for  which  solutions  are  known.  The  impedance  Z  consists 
of  three  terms:  (i)  the  radiation  resistance  R  due  to  the  propagating  TEM 
mode,  (ii)  the  reactive  ordinary  rectangular  waveguide  contribution  Z-^,  and 
(iii)  the  gap  field  perturbation  reactance  principal  value  integral.  Ideally, 
since  each  of  these  terms  has  a  distinct  physical  interpretation,  each  term 
should  be  checked  individually.  Because  the  probe  is  electrically  short,  the 
probe  radiation  resistance  may  be  compared  to  that  of  a  Hertzian  dipole  via  a 
result  given  earlier  by  Tippet  [1].  His  expression  for  the  radiation 
resistance  of  a  Hertzian  dipole  involves  the  square  of  IEQdJ£  where  I  is  the 
current  magnitude,  EQ  is  the  TEM  mode  electric  field  distribution,  and  di   is 
the  effective  dipole  length.  In  the  monopole  case,  the  dipole  quantity 
(IEQdJi),  may  be  replaced  by  an  effective  value  which  results  from  integrating 
the  product  I(y)  EQ(y)  over  the  monople  length.  The  two  results  agree  up  to 
the  order  of  our  approximations.  Next,  by  letting  the  gap  width  go  to  zero 
(g+0)  thereby  reducing  the  upper  chamber  to  an  ordinary  rectangular  waveguide, 
Z-l  may  be  compared  to  a  result  given  by  Collin  [3].  If  we  allow  the  frequency 
to  be  large  enough  that  the  ordinary  rectangular  waveguide  TE^q  mode  will 
propagate,  then  our  result  reduces  exactly  to  that  found  by  Collin.  Both 
these  checks  are  detailed  in  Appendix  D  for  the  special  case  when  a^=a2=a. 

The  perturbation  reactance  term  is  more  difficult  to  verify  and  is  not 
comparable  to  any  simple  special  case  result.  However,  a  good  check  of  the 
self  consistency  of  our  expression  exists.  If  we  had  allowed  k  to  be 
sufficiently  large,  such  that  k2>M^,  the  expression  (12)  for  Z^  would  contain 
poles  at  cx=+bIq  ,  where  &!Q  =  (k  -M^r2  .  This  would  imply  an  apparent  mode 
propagating  the  same  as  an  upper  chamber  TE^q  mode.  Equation  (15)  for  Z2 
would  also  imply  the  existence  of  this  mode.  However,  the  guide  as  a  whole 
supports  no  such  mode,  therefore,  we  expect  that  the  contributions  from  Z^  and 
Z2,  to  this  evident  Te[q  mode,  should  cancel  out  to  the  order  of  our  previous 


approximations.  It  is  found  in  Appendix  D  that  the  radiation  resistance  due 
to  a  propagating  TEJq'  mode  in  a  rectangular  waveguide  of  the  dimensions  of 
the  upper  chamber,  is  given  by 

n     no     .  2  kd  ,„, 

R  =  pry  tan  —7y  (23) 

2»ibikBiS 

Now  consider  the  perturbation  contribution  given  by  Z2  in  (15).  In  order  to 

avoid  the  degeneracy  due  to  an  additional  TEJq'  mode  also  propagating  in  the 

lower  chamber,  we  refrain  from  letting  the  chamber  widths  be  equal  until  after 

we  have  evaluated  the  pole  contribution  at  a  =  +$!n  .   As  a  approaches 
l\\  (\\  ( i  \o  o    1  /2       iU 

tPiJ  ,  we  see  that  V.    '=   (b{q   -  a  )  '     +0.  Evaluating  the  half  residues 

at  a  =  ±$10  ,  we  find 

2a1b1(31(n)2  "  <?) 
L(0)  +  l   l     1U (24) 


£   +  (7i/2a1)2  (25) 


2aibik  (3lo  -  a)  l 


JoEtfM^+a2)1^  +  JQ(kt)  (27) 


Because  both  irg/2a1  and  kt  are  small,  the  Bessel  functions  may  be  replaced  by 
unity.  Collecting  results  gives 


h 


-lwyk       csc  kd  -2th)      in  1)  7— zr — ry  h  ,  (e._   ,d) 
(1)--— 2  «*±31Q      U72i^       1     10 


a=±31Q  4ax 


(28) 


"no  .2  kd 


2alblk  310 


or    tan  ~2 
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which  indeed  will  cancel  the  unperturbed  result  given  by  expression  (23). 
Thus  the  poles  at  a  =  ±3J0  ,  that  apparently  imply  a  physically  meaningless 
mode,  do  not  contribute  up  to  the  order  of  our  approximations.  The  impedance 
expressions  will  next  be  used  to  generate  numerical  data  for  typical  TEM  cell 
geometries. 

The  radiation  resistance  R  given  by  (19)  is  plotted  in  figure  3  for 
g/a=0.1  and  0.2.  The  dimensions  of  the  RCTL  cross  section  are  a^=a2=1.0  m, 
and  b-j_=b2=1.0  m.  The  probe  varies  in  length  from  d/a=0.0  to  d/a=0.95,  and  is 
of  radius  t=l  mm.  The  frequency  is  set  at  1  MHz.  The  curves  show  a  basically 
parabolic  dependence  on  d  with  the  main  difference  between  the  two  curves 
resulting  from  the  logarithmic  dependence  of  Zc  on  the  gap  size.  This  d^ 
variation  is  significantly  larger  than  the  free-space  probe  radiation 
resistance  which  is  proportional  to  (kd)  ,  a  frequency  dependent  quantity. 
This  difference  between  the  radiated  power  by  a  short  monopole  in  free  space 
and  that  delivered  to  a  TEM  mode  has  also  been  encountered  in  other  types  of 
problems  such  as  the  coupling  between  a  dipole  antenna  and  an  infinite  cable 
[8].  Figure  4  gives  the  dependence  on  d  of  the  ordinary  rectangular  waveguide 
reactance  X.  The  same  set  of  parameters  is  used  to  plot  X,  given  by  (20),  as 
that  used  to  plot  R.  The  curve  shows  a  highly  capacitive  contribution  which 
is  independent  of  the  gap  size.  Again,  this  result  may  be  compared  to  the 
free-space  result  for  an  electrically  short  monopole  which  is  known  as 
proportional  to  -i(kd)"1  [9].  From  (20),  the  dominant  reactance  behavior  for 
small  kd  is  more  like  -i(kd)  .  Thus,  we  find  that  the  reactance  of  a 
monopole  in  a  TEM  cell  is  more  highly  capacitive  than  that  in  free  space. 
Finally,  the  perturbation  reactance  6X  is  plotted  in  figure  5  using  (21), 
again  for  the  same  parameter  set.  We  find  that  |6X|<lfi  for  the  entire  range 
of  probe  lengths.  Thus  6X  is  at  least  three  orders  of  magnitude  less  than 
X.  Therefore,  neglecting  6X  will  yield  excellent  results.  The  FORTRAN 
program  used  to  compute  Z  is  listed  in  Appendix  E. 
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6.  CONCLUDING  REMARKS 

This  report  has  detailed  the  manner  in  which  a  coaxial  probe  excites  a  TEM 
cell.  The  analysis  utilizes  three  simplifying  assumptions,  which  are  (i)  the 
gap  is  small,  that  is  G  £nG<<l  where  G=irg/2a,  (ii)  the  probe  is  thin, 
or  irt/2a<<l  where  t  is  the  probe  radius,  and  (iii)  the  frequency  is  not  too 
high  so  that  only  the  TEM  mode  propagates  and  (ka)^  is  small  when  compared  to 
unity.  None  of  the  requirements  is  restrictive  in  terms  of  common  TEM  cell 
usage.  The  resulting  expression  for  the  probe  input  impedance  consists  of 
three  distinct  terms  R,  X,  and  6X,  these  being  the  radiation  resistance  due  to 
the  propagating  TEM  mode,  the  ordinary  rectangular  waveguide  reactance,  and 
the  reactance  perturbation  due  to  the  gap  fields.  The  first  two  terms,  R  and 
X,  are  computationally  simple,  with  R  showing  a  basically  <r   dependence  on 
probe  length  and  X  giving  a  highly  capacitive  reactance.  Although  5X  consists 
of  a  nontrivial  principal  value  integral,  numerical  results  show 
that  |6X|<<X,  and  therefore  6X  may  be  neglected.  Thus  a  reasonably 
straightforward  form  for  Z  is  achieved  which  may  be  used  to  analyze  probes 
inserted  into  a  TEM  cell  either  to  measure,  or  excite  fields  as  well  as  other 
types  of  feed  lines. 
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APPENDIX  A 

DERIVATION  OF  k{^  and  K^ 

Generalizing  the  chamber  widths  involves  some  modifications  in  the 
expressions  derived  in  NBS  Technical  Note  103  [4].  For  the  most  part,  this 
involves  nothing  more  than  replacing  a  by  a..  The  details  will  be  briefly 
outlined  beginning  with  the  TE  type  mode  kernel  which  becomes 

Q(h)(t,f )  =  Q^h)(t,f)  -  Q|h)(t,f)  (a-1) 

Q^h)(t,t')  s  -I*n  [(^(4^)  It2  -  t'H]  (a-2) 


Q^h)(t,t')  ,  f  1  I 
j=l   J   mo 


cot  rJ'V     -, 
JO J_  +  J_ 


K(j)         Mj 
m  J 


=   Qf        (a-3) 


where  t  =  a-x,  t'  =  a-x',  a  =  min(al5ao),  V^ '  =  c2  -  M.»  M.  =  mTr/2a., 
and  mo  refers  to  the  summation  over  odd  m.  This  kernel  will  result  in  the 
following  solution  for  the  unknown  gap  field 

f(h)(f)  =  -Id* L(a)f!(x')  (a-4) 

2a1/g2-t'2       l 

^-  M^J  +  h]"1  (a-5) 

M.  cos  K^(b.  -y) 

f.(xj  =  I         -1  mm    fa         sin  ^  cos  Mx  J  (Mg)       (a-6) 

J  z         mo     Ku;sin  KUV      ^     J   °  J 
m      m   j 

where  we  recall  that  x!  is  the  dipole  location.  The  above  solution  assumes 
that  the  dipole  is  centrally  located. 

The  y-component  of  the  electric  field  in  the  upper,  and  lower  chambers 
may  now  be  found.  Both  TE  and  TM  type  unperturbed  (ordinary)  modes,  as  well 
as  TE  type  perturbed  modes  will  contribute  significantly  to  E  J'  .  The 
unperturbed  modes  exist  only  in  the  upper  chamber.  Denoting  the  unperturbed 
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contribution  by  RiJ '    ,   and  the  gap  perturbation  contribution  by  Ko  ,  and 
referring  to  Appendix  C  of  NBS  Technical  Note  1037  [4],  we  find 


Kl   Ut'V   6ij   ^  3x  z  Ut'  V   c2  3y  Ez  Ut'V 
Suitably  modifying  equation  (49)  of  NBS  Technical  Note  1037  [4]  yields 


(a-7) 


K^^,  Xt)       ^   ^ 


2aiaC 


2~  Idfi,  L(a)  f  (xt)  f1  (xt) 


U-8) 


Taking  the  proper  derivatives  required  in  (a-7)  and  substituting  the 
proper  forms  for  E^  ,  and  H^ '     allows  K,^'  to  be  written  as  follows. 

AnMj  +  4  N  j 
K|j)(xt,x;)  =  6   M  Id£   J    \     k   3   g(j)(xtJx;)        (a-9) 
1    '  r     U  C     m,n=0   M.  +  N^  -£   mn     r 

■J       U 


g   (x.,x')  =  ( — r— )  sin  M.(x  +  a,-)  sin(x'+a.)  cos  N  .y  cos  N.y'    (a-10) 
mn   t  V       va_.b_.      v    j'        v     3      J 


J  J 


where 


A  = 
n 


\     n  =  0 


1  n  >  1 


(a-11) 


and  N.  =  mr/b..  It  will  prove  convenient  to  extract  the  delta  function 

J   ~(i)J  -   -• 

nature  in  IOJ/  which  occurs  when  x.  ->  x  . 


lotice  that 


A  m2  +  -2L  N2  =  A  \  [M2  +  N2  -  ?2  +  %  M2  +  C2] 
n  J   7  J    n  k2   j    j       a2  j 


(a-12) 


Therefore  K;   becomes 


Kij)(x+,x:) 


1  -t^-lj^^^^/n^^) 


I   9  J  ■>      ?  g(j)(x.,x:) 

m.n=0  M2^  -  s2  """>   l  l 


(a-13) 


If  we  expand  a  delta  function  in  a  double  Fourier  series  over  a  chamber  cross 
section  we  will  find 


m,n=0 


(a-14) 


15 


Therefore  K 


(J) 


1 
(j) 


becomes 


icuy 


Kl   Ut'  V  "  6lj  ,2 


Id£ 


a 


C 


,  «(xt 


*;> 


°°   A  [M2+a2]  /  -v 

X    "  Jo   o   o  9(j)(x+,  x  ) 

m,n=0    M2+N2-c2  mn   l   l 
J   J 


(a-15) 


R{J> 


and  K 


(j) 


represent  the  y-component  of  the  electric  field  excited 


by  a  symmetrically  located  VED  at  x,  .  There  remain  only  two  small  changes 
to  adapt  this  result  to  the  probe  excitation  problem.  First,  the  probe 
carries  a  negative  y-directed  current.  Thus  Id£  -*  -  idl   .  Second,  the 
Green's  function  represents  results  from  a  unit  strength  dipole.  Therefore 
-Id£  ->  -1 


Therefore  K,   and  Ko    finally  become 


i  'vV " 


lj  k2 


A  [M2  +  a2] 
nL  j    J 


m,n=0  M2  +  N2 


C 


V(*t-xP 


g(j)  (x.,  x!) 
3mn    t   t 


K2  (xt'  Xt} 


lojyiT 


2al3jC 


2  L(a)fJ.(xt)f1(xt) 


(a-16) 


(a-17) 
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APPENDIX  B 


EVALUATION  OF  Z 


The  integral  defined  for  Z,  is  to  be  evaluated. 


zl  =772 


2ttI 


'V\    i  uO)iv    v<\J°<z 


t    s  s 


J(x)       K)u(xt,x^)eiazJ(x,)dads'  ds  (b-1) 


(1) 


Kr;(x,,  xj  =  -  Jwy 


,2  ?      Ut     Xt    j 


2       2 

Mi +a        rn  -  -i 

m;n=0       n     M2+N2-?2     mn       l     l 

The  a-integration  may  be  performed  first,  which  requires  evaluating  the 
following  pair  of  integrals. 


(b-2) 


I 


I, 


2 
a  .,—    — . ,  iaz, 
-~2  o(x.  -  x  )e   da 

2   2 

I    A    1     .  .  g(1)(xt,x:)elaz  da 
m,n=0  n  M2+N2-C2   mn   t  l 


(b-3) 


(b-4) 


In  each  integral  we  will  deform  into  the  upper  a  half  plane  for  z  >  0 


(forward  waves),  and  into  the  lower  a  half  plane  for  z  <  0  (backward  waves). 

iklzl 


For  I, ,  we  find 


I,  =  -iirk6(x.f  -  x.  )e 

The  denominator  of  I~  may  be  factored  as  follows 

m2       m2       r2         2       r(l)2 
11  mn 

T(1)   =   [M2  +  N?  -   k2]2 
mn         L   1  J 


(b-5) 


(b-6) 
(b-7) 
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Thus  we  see  that  a  simple  pole  arises  in  the  I_  integrand  when  ct  -  +ir 

2    3  mn 

Again  deforming  as  for  I,  gives  the  following  result. 

k2-N?  m       -r(1)|z| 

l2  =  *       I       ^-TTT^Vxpe  ^ 
m,n=0  n  ru;   mn   L  z 
mn 

We  next  proceed  to  the  integration  over  the  current  distributions 
J(x)  and  J(x').  The  source  current  J(x')  will  be  concentrated  on  an 
idealized  linear  antenna  located  at  the  origin,  i.e.  x'  =  (0,y ' ) .  There- 
to 


(b-8) 


re  x.  f   x, i  on  S  and  the  delta  function  term  will  not  contribute  to 


Z-,  •  Thus  Z-,  reduces  to 


I     =  ML 


o,272  _  L_n   "n  JJT 


1  2  M2 
k  -  N-, 


2k  I"  m,n=0 


r 


J(x)g(1)(x,,x:)j(x')e 
v  /3mn   t  V    v  ' 


mn  ' 


mn 


S  S 


ds'ds 
(b-9) 


The  simpler  integration  over  the  source  current  will  be  considered  first, 

We  recall  that  the  assumed  current  density  on  the  probe  will  be 

I 
J(x)  =  ^t  sink(b!  -d-y)6(p-t) 

For  the  6*  integration  we  let  t'  +  0,  and  therefore 


(b-10) 


-r(1)|z 
I3(^-)I0  sin  M^x  +  a^  sin  f-  cos  Nny  e 


'V 


(b-11) 


where  the  integral  I~  is  given  by 

Io 


fbl 


brd 


, ,       ,      ,  x        ..     ,,  .  ,        k  cos  n-n-(cos  kd  -  cos  f^d  )    ,,    19\ 
sin  k  (b.-d-y'jcos  N,  d' dy'    =  ^ ~ 1 — *-  (b-12) 

1         '  k  -N^ 


Because  of  the  factor  sin  rmr/2,  only  the  summation  over  odd m  will  contribute. 
•,  ta 

loop  I 


Therefore  Z-,  takes  the  form 


A„  cos  niT 


Z,   =  — = %-        I        — — pry-  (cos  kd  -  cos  N,d)      J(x)cosM,xcosN,y 

alblkIt      m0'n=0         rin  M;  JS 


-r(1)|z| 
xe     mn         ds 


(b-13) 
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The  integration  over  S  may  be  split  into  separate  integrations  over  y 
and  0  .  The  y-integration  has  already  been  encountered  and  we  find 

I  Io  P  -riPtlsine 


-r^iz 


J(x)cos  M,x  cos  N-,y  e 


mn   '   ds  =  ~-%-4-  I     cos(M,t  cose)e  ~mn  t  de 

2fTt   J      v  1        ' 

o  (b-14) 


If  we  now  interchange  the  remaining  theta  integration  with  the  double 
summation,  Z,  may  be  written 


zi = 


ioou  I' 


Tra1b1It 


(cos  kd  -  cos  N-j d) 


n=0 


k2-N^ 


x  I        ,,  x  cos(M-.t  cos  6  )e 


■r^t  sin 
mn 


mo  r 


mn 


de 


(b-15) 


where  we  have  noted  that  the  theta  integration  is  symmetric  about  tt  .  Let 
us  now  consider  the  sum  over  m 

-r^  'tsin  6 


S  (9)  =  I        pjy-COs(M,tCOSe)e 

mo  rmn 


mn 


(b-16) 


The  dominant  series  will  be  for  n  =  0.  S  (e)  may  be  expressed  as  follows. 


S  (0)  =  Re  I 
mo 


For  large     m     we  see  tha 


iM,t  cose  -  rijt  sini 
e  e 


r 


nr 


mn 


(b-17) 


n=0 


t    r'1'- 


'^   M, .  Adding  and  subtracting  the  asymptotic 
mo    1       3  J 

form  of  the  summand  will  thus  give 


S  (0)  =  Re  I 
mo 


+  Re  I    e 
mo 


M,t(i  cos0  -  sin  0) 


iM,t  cos 


-r^t  sin0 


mn 


■M,t  sin  0 


mn 


'1   Jn=0 


(b-18) 


where  the  first  series  represents  the  major  contribution.  The  first  series 
may  be  summed 
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I 

mo 


iM^e19) 


•a,        ,  10     ia, 
-Un  [tan(^f-)]  +  ^i 


4a 


(b-19) 


We  now  assume  that  the  probe  is  thin,  irt/4a,  <<  1,  which  allows  us  to 
replace  the  tangent  by  its  argument 

i  ft 

an[tan(^|f-)]z  ^   (^  eie)  =  in   (jjj-)  +  ie 
If  we  now  take  the  real  part  of  the  first  sum  in  (b-18)  we  see  that 


(b-20) 


a,    4a,         iM,tcos0 
Sn(6)  ~  ln{-±)   +  Re  I    e 

mo 


TTt 


(1). 


-Tv  't  sin  0   -M,  t  sin 

«  mn        ~  1 
e  e 


r 


nr 


mn 


The  correction  series  will  be  considered  next.  Again  consider  F 


"1 

(1) 
mo 


Jn=0 
(b-21) 


rtio)=  Ml(1  "k2/Ml)2  *   V1  -k2/2^)  =  V1  -£)  (b-22) 

2   2 
where  e  =  k  /2M,   is  assumed  to  be  small.  The  correction  series  may  be 


written 


i 

mo 


M,  t(i  cos0  -  sin  0  ) 


1 


M,  t esin  0 

ej -  1 

1  -e 


(b-23) 


■1 


For  e  small,  we  have  (1  -  e)   «  1+e.  This  series  converges  very   rapidly, 
thus  only  a  few  significant  terms  may  be  considered.  Because  Trt/2a,  «1 
for  these  few  terms,  the  exponents  may  be  replaced  by  unity  with  little 
error  giving 

,  2 


y  —  -   K     y 

^  M,  ~  T  L      M3 
mo  1     c     mo  M, 


a   3,2 
4a,  k 


I 


} 


(b-24) 


mo  m 


The  sum  of  the  inverse  cube  over  odd  m  may  be  done  in  terms  of  the  Riemann 
Zeta  function.  Recalling  that 


C(s)  =  I     k 
k=l 


-s 


(b-25) 


For  Re(s)  >  1 ,  we  see  that  our  expression  may  be  written 


l 


1 


I 


1 


mo  m 


-  8  c(3) 


(b-26) 


m=l  m    m=l  (2m) 
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The  Riemann-Zeta  function  is  well  tabulated,  and  from  [6]  we  find 


c(3)  z   1.20205 


(b-27) 


Collecting  these  results  we  have  the  following  for  the  dominant  series 

s0(e) 


3  2 
a,     4a    a,k 

S  (6)  =  -L  In    (-T1)  +  -^5-  (4.207175) 

O         TT        TTt  3    V  ' 


(b-28) 


The  remaining  series  S  (0)  for  n  >  1  will  be  summed  in  a  manner 
3        n        — 

similar  to  that  used  by  Collin  [3]  in  treating  the  probe  excitation  of  a 
rectangular  waveguide.  We  will  make  use  of  the  following  Poisson  summation 
formula:  if  f(x)  is  a  square  integrable  function,  then 


I     f(na)  =  -         I     F(^) 
n=  -o°         n=  -°° 


(b-29) 


where     F(co)      is  the  Fourier-transform  of     f(x)     defined  by 


IcoX 


F(co)   =       f(x)e       dx 


(b-30) 


Our  particular  problem  involves  the  function  f(x)  defined  by 

ixt  cos  9  -1^  (x)t  sin 


f(x)  =  e 


mn 


/r(1)(x) 
mn  v  ' 


T(1)(x)  =  [x2  +  N2  -  k2]2 

mn   V         L 


(b-31) 
(b-32) 


This  function  is  clearly  Lp  so  Poisson's  summation  formula  applies.  The 
Fourier  transform  is  given  by 


F(o>) 


0°  -r^'(x)t  sin  e 

'  e  mn         ix(tcose+o)) 

»     mn 


(b-33) 


,(D 


If  we  now  rewrite  the  complex  exponential,  and  notice  that  r^  (x)  is 


even  in  x,  then  F(w)  takes  the  form 


°o  -r^  (x)t  sine 


F(uO  =  2 


fe  mn 


mn 


cos((t  cos  0  +  co)x)dx 


(b-34) 
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This  integral  may  be  evaluated  in  terms  of  the  modified  Bessel  function 

K  [7]  with  the  result 
o 

7  V 

F(co)   =  2KQ[(N2-kV   (t2sin29  +   (tcose  +  w)V]  (b-35) 

If  we  now  return  to  the  definition  of     S   (0)  we  see  that 

n 

s(e)  =  1  +l       f(^)  =  1  I     If(^L)  -  f(f )}  (b-36) 

Thus  applying  Poisson's  summation  formula  we  have 

S(e)=|   I   {-J-F(4a1m)--J-F(2a  m)}  (b-37) 

m=  -oo 

which  in  terms  of  modified  Bessel  functions  K   becomes 

o 


k  o  ok 


S  (0)  =  —        J       2K  [(N?  -  K2)^((4a1m)2  +  8anmtcos  e  +  t2)2] 
n  ir  ^     j     ol  1  1 

m=  -°° 

-  Kq[(N2  -k2)**   ((2a1m)2  +  4a1mtcos  e  +  t2)^]  j  (b-38) 

The  modified  Bessel   function  decays  very  rapidly  so  the  only  significant 
contribution  comes  from  m=0.     Thus 

Sn(9)   ,^Ko[(N2-k2)\]  (n>l)  ^  (b-39) 

Returning  to  the  integral   for     Z,    in   (b-15)  we  have  the  following 

2         TT  2 

icou  I       f  (/,  .  ,n2  °°     (cos  kd  -  cos  Nnd)  \ 

V— H         (         2 Skd)      So(6)+     I     ' 2 r1 Ve>de 

1   ]   t    ol  1  J       (b_4Q) 

Note  that  the  approximate  form  for     S   (0)   in   (b-39)    is  now  independent  of 
0.     Thus  combining  all   these  results,  we  have 

2  /  2  2 

icoy  r  0  4a,  a^k" 

Z,  = ^V  (1  -  cos  kdT      In   {—±  )   +  -~-  (4.207175) 

'        Zvb^l[  |  TTt  ^2 

-     -     /       sin2N1d/2  x2      K  [t(N2-k2)%] 

-2k2     I       1 ^l—  )     — i1-^—  (b"41) 

n=l   V       sin  kd/2    /  N^  -  \C 
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Note  that  Z,  given  in  (b-41)  is  approximate,  based  on  two  assumptions 
that  the  probe  is  thin  (Trt/2a-,  «  1)  and  that  the  frequency  is  not  too  high 
(kaj  «  1). 

Finally,  we  note  that  the  total  current  at  the  aperture  (y =  b, )  is 

given  by  I.  =  -  I  sinkd.   Thus, 
to 

2  2 

1     z  ^M     tan2   Win(i)+  i>      (4.207175) 

1        2TTb,r  l  irt  /  (b_42) 


1 


o  /  sin2  N,    d/2   \         K  [t(N2  -  k2ft 

2k       I     V-— "H ^ 12 2- 

n=l    ^  sin     kd/2       '  N^  -   VT 
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APPENDIX  C 


EVALUATION  OF  Z, 


The  integral  defined  for  1?     in  (13)  is  to  be  evaluated 


Z2  "  „_T2 


2ttI 


0)  ^   IT'  \oiaz 


J(x)   lO,u  (xt,  xt  )e'uzJ(x  )dads'ds         (c-1) 


t  S  S 


where 


and 


41)(\'  \]    "TTZ  ^^V^V't1 


2a1? 

'8/a,a0\     2 


frg 


a  . 


j  =  l  j  mo  L   K 


cot  K^b.     , 
(j)       M, 


W  =   I   (1)  "  0  .    sln^costyj^g) 

mo  Kv  ;sin  V    'b, 
m      m   1 


(c-2) 


(c-3) 


(c-4) 


If  we  define  two  integrals  I.  and  l^     via 


— i  \  i  i 


J(x  )  f1(xt)ds 


(c-5) 


J(x)   f1(xt)eiotZds 


(c-6) 


and  interchange  the  a-  and  surface  integrations,  we  may  write  Z~  in  the 
following  form 


7    -lojy 
?  2  2 

4al!t 


'  Uo) 


2  !4  h   da 


(c-7) 


The  integrals  I.  and  I,-  now  need  to  be  evaluated, 
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The  assumed  current  distribution  is  given  by  (2).  Therefore  I-  becomes 
V1"  TJo(Ml9)  ^ 


h  =  I       I 


'o    ±        M)ein  „(U 


mo         10    ysin  IO    yb, 

m  m       1    b, -d 


sink(b1  -d-y')   cos  K^1)(b1  -y')dy 


(c-8) 


Performing  the  y-integration  gives 

rmr 


4  o 


M 


1  sin-^-  J0(M-,g) 


K^   sin  K^   b, 
m  m 


(cos  kd  -  cosK^d) 

m 

M2    X       2 

M,    +  a 


(c-9) 


The  integral   Ir  is  similar  to  I,  except  that  an  additional  azimuthal 
integration  around  a  probe  must  be  performed.  The  y-integration  is  the 
same  as  in  I,  so  Ir  becomes 


kl 


I 


-  y 

rt  L 


M.sin  ^  J   (M,g)  (cos  kd  -  cos  K^d) 

1  2      oM3/  v  m       ' 


mo       Kv    'sin   10    'b,  M,   +  a 

m  m       1 


x    cos(M1tcose)e1atsin9tde  (c-10) 

o 
If  we  designate  the  theta  integration  by  H  (t),  then  letting  0  -»-  2tt  -  9 

in  the  interval  [tt ,2tt]  allows  H  (t)  to  be  written 

H  (t)  =  2t   cos(M,t  cos  9 )cos(at  sin  6  )d9  (c-11) 

o 
The  second  cosine  term  may  be  expanded  in  terms  of  Bessel  functions  [7] 

oo 

cos(at  sin  e)  =  J.(at)  +  2  \     J9. (at)cos  2ke  (c-12) 

0        k=l  ^ 

Interchanging  the  sum  and  integral  gives 

H  (t)  =  2t  jj  (at)   cos(M]t  cos0)de 

0         TT 

f  ) 

+  2     I     J9,(at)       cos  2k6cos(M,t  cos  e)dej  (c-13) 

k=l  J 

Again  consulting  [7]  allows  us  to  evaluate  these  integrals  with  the  result 
that 

Hm(t)  =  2Trt{  Jo(at)JQ(M1t)  +  2  I      (-l)kJ2k(at)J2k(M1t)J      (c-14) 
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This  is  of  the  form  described  by  Gegenbauer's  Addition  Theorem  [6],  which 
for  the  special  case  that  the  function  being  expanded  is  a  zero  th  order 
Bessel  function,  becomes 


J    (pR)   =  J    (pq)J    (pr)   +     I     cos   k<j>  J.  (pq)J.  (pr) 


o 


o 


(c-15) 


k=l 


7  7  7 

where     R  =   [q     +  r     -   2qr  cos  $  ] 2.      If  we  let     p  =  t,  q  =  a,   r  =  M, ,   and 

T7 


<p    =  j  ,  we  f  ind 


R  =   [a2  +  M2]^ 


(c-16) 


Therefore,  H  (t)  in  (c-14)  becomes 

m       v    ' 

Hm(t)  =  2^tJQ[t(a2  +  M2)%] 


(c-17) 


If  for  simplicity,  we  define  the  following  quantity 


h  (a,d) 


M1  sin^  J0(M19)    (coskd  "cos  Kmd}) 


rnir 


V**   sin  K^  b, 

m      m 


2    2 

1^  +  a 


(c-18) 


the  above  results  for  I.  and  Ic  combine  to  give  Z0  as 

4       5  3     2 


h  = 


2  2 

-iwyk  I 
0 

a   2  T2 
4al  Jt 


^   T   h  ,(a,d)  I    h  (a,d)J  [t(a2+M2)%](d 
2     _  •  m       „.,  m     o      1   1 


mo 


mo 


a 
(6-19) 


If  we  now  recall  that  I  ,  the  total  current  at  the  aperture,  is  given  by 

-I  sin  kd,  Z0  becomes 
o        2 


7   -iojyk  „„Jl.   ■  [  \   L(ot) 
Z9  = £ —  esc  kd 

4ai 


2.„2^- 


f  I  hm,(a,d)  £  h  (a,d)J  [t(aSM^)  ] 
l,        mo        mo 


da 
'(c-20) 


Unfortunately,  the  a-integration  is  sufficiently  complex  that  some  numerical 
approach  has  to  be  used  in  evaluating  this  integral.  Notice  that  for  large 
m,  since  d  <  b-, ,  h  (a,d)  will  contain  an  exponential  decay  according  to 
exp(-a(b-,-d)).  Therefore,  the  integrand  will  decay  rapidly  with  a  and  a 
numerical  routine  should  not  prove  time  consuming. 
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APPENDIX  D 
CHECKS  OF  R  AND  Z 


This  appendix  will  detail  checks  of  two  of  the  impedance  contributions. 

First,  the  radiation  resistance  will  be  checked  against  a  result  given  by 

Tippet  [1],  and  second,  Z,  will  be  compared  to  the  rectangular  waveguide 

result  given  by  Collin  [3].  The  radiation  resistance  r'  of  an  electric 

dipole  inside  a  TEM  cell  due  to  the  TEM  mode,  as  given  by  Tippet  [1],  is 
2Z      IE  cose  2 
Re  =  — \  (d£   °2V  }  (d_1) 

where  Z   is  the  characteristic  impedance  of  an  RCTL,  I  is  the  magnitude 
of  the  current,  d£  is  the  effective  length  of  the  dipole,  E   is  the  TEM 
mode  electric  field  distribution,  9  is  the  angle  between  the  dipole  direction 
and  the  y  axis,  (for  a  vertical  dipole  9=  0),  and  V  is  the  voltage  asso- 
ciated with  the  TEM  mode.  Because  the  probe  is  an  electrically  short  mono- 
pole  (kd«l)  the  radiation  resistance  should  be  nearly  that  of  a  dipole 
with  the  same  current  distribution.  Expression  (d-1 )  assumes  an  elementary 
dipole,  however.  We  may  replace  IE^d£  with  an  effective  value  defined  by 

(IE  d£)  ,, 
v  o   eff 


o 

h 

I(y)E(y)dy  (d-2) 

-d 


where  I(y)  =  sin  k  (b  -d-y),  and  E(y),  for  a  centrally  located  dipole,  is 
found  in  [4]  to  be 

?u     cosh  M(b,-y) 

E(y)  =  —  I     -tJ sin  ^  J  (Mg)  (d-3) 

VJ;    a   L         sin  h  Mb,       2   ov  a/  x   ' 

mo        1 

Integrating  yields  the  following 

(IE  d£)  «  =  2B.  i    (cos  kd-  cosh  Md)  sin  rmr  j  (  }       (d_4) 
v  o  yeff    a   L  sin  h  Mb,         2   ov  3/ 

which,  if  we  recall  the  definition  of  hm(a,d),  may  be  written 


m 
mo 


(IE  da)    .-   =  -^2Vk.  y  h  ((<  d)  (d_5) 

v  o  7eff    a  Jtl    mx  ' 
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If  we  now  normalize  by  the  magnitude  of  the  current,  |l|  =  sin  kd,  we  find 


R   becomes 

2Z  k' 

K  =   — >r  csc^kd  (I    h _(k,d)f  (d-6) 

a  mo  m 


97  k2   9  2 

c     2,  -  ^ 


This  may  be  compared  to  the  radiation  resistance  R  due  to  the  TEM  mode 

given  by 

2Z  k2 

R  =  -f-csc2kd  £  h  ,  (k,d)  I     h(k,d)J  [t(k2+M2)^]  (d-7) 

a         mo  m      mo  m     ° 

The  term  h  (k,d)  has  a  strong  decay  in  M  according  to 

hm(k,d)^e    '  (d-8) 

M 

Only  a  few  terms  will  be  significant.  For  these  few  terms  J   will  remain 

J  3  0 

near  unity  since  we  have  assumed  that  7Tt/2a  «  1.  Thus,  the  radiation 
resistance  may  be  approximated  by 

U/         ?      I  \2 

R  =  — V-  csc^kd  I    h  (k,d)  =   R'  (d-9) 

a        \mo       / 
and  the  results  indeed  do  check. 

Next  the  case  of  no  gap  will  be  considered.  As  g  -t  0,  L(a)  ->  0 
so  only  Z-,  contributes.  Collin  discusses  a  rectangular  waveguide  with 
a  short  circuit  at  a  distance  i     from  the  probe  [3].  Thus  we  need  to 
modify  his  result  by  letting  £  ->  °°.  Making  the  appropriate  changes,  the 
radiation  resistance  for  a  probe  in  a  rectangular  waveguide  becomes 

R  ■  2ib\k  ta"2  f  <d-10> 

The  equivalent  reactance  is  given  by 

X  =  2^kta"2¥  U(£)+^(.2072)-2(l-i) 

'  ^  k 

7  I         sin2N1d/2?Kn[t(N2-k2)2]         ) 

-2k2  I     fi-_l    )-5 1  (d-11) 

n=l  \    sin  kd/2/      N^  -  r    ) 

In  order  to  compare  our  expression  to  that  of  Collin,  we  must  return  to 

Appendix  B  (b-4)  and  let  k  be  larger  than  M  (a-.  =a2  =a).  This  will  create 
a  pole  on  the  real  axis  when  a  =  3-in-  Thus  I?  becomes 
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2         2B1Q  910  lV     t'e 


^2-n?  -r<1>Ui 

m,n=0    "     TKU       "in     t     x 

mn  (d-12) 


where  the  primed  summation  indicates  that  the  term  m  =  "I ,  n  =  0  is  missing, 
Performing  the  current  integrations  yields 
.2 


Zl   = 


o 

nab  I2 

t      o 


(1  -  cos  kd)    „„   /Tit  o       Av 
—^ —L-  cos(~—  cos  9)e 


i 3  Qt:  sin 


2ik2e 


10 


°°  (cos  kd  -  cos  Nnd)' 

t      A    - 

^  n 


mo,n=0 


2       .,2 


k     -  N 


-7-yyC0s(Mt  cos  0)e 


-ril't  sine 


mn 


de 


mn 


(d-13) 


The  sum  over     m     and     n     is  precisely  that  encountered  by  Collin   [3]  with 
the  result  that 

,  -r^'tsine 


I  A 

L.         _        r 


mo,n=0 


(cos  kd  -  cos  N,d) 
~~     2       2~ 


yy  cos(Mt  cos  e)e 


mn 


mn 


2 
_   (1  -  cos  kd)        a_ 

2k2 


2  2 

in(™)  +  i<_A_  (.2072)  -2(1    -  -) 
Tit  2  a 


v     (cos  kd  -  cos  N-,d)' 
+  i.      L     


n=l 


k2-N2 


KQ[t(Nf  -  kZ)    ] 


The  significant  contribution  for     Trt/2a  «  1     may  be  approximated  via 

1 


,  i3int  sin 

cos(^  cose   )e     IU 


(d-14) 


(d-15) 


Combining  these  results  and   integrating  over  0   ,     which  simply  introduces  a 
factor  of  tt,  yields 
.2 


Zl    - 


iwy  I( 
-nab  I! 


(1-coskd)2^  a(l-coskd)2'      rn(4i}  +  kY(t2072) 
2ik  3,~  2k^  / 


■2(1  --) 


10 

2  P 

<*>  (cos  kd  -  cos  N,d)      .,  r.,   2     ,2N2-, 

+  a     I ,       ?]         Ko[t(N  1  "  k  }   ] 

n=l  k    -N^ 


(d-16) 
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Finally,   letting     I     =  -I     sin  kd,    and  separating     Z,      into  real   and  imag- 
inary parts     R     and     X     respectively,  we  obtain 
n, 


R  = 


X  = 


o 


2abk3 


n, 


10 


.      2  kd 
tan    -=- 


'o       .      2  kd 
^Tbk     tan     ~1 


2  2 

£n(%  +  M-  (.2072)   -   2(1  -  ^) 
v7Tt  2  a 

7T 


(d-17) 


_2k2    ?L     sm^d/zx2     KQ[t(N2-k2)2]       j 
n=l\       sin2kd/2    /  N2  -  k2  J 

which  are  the  precise  results  found  by  Collin  in   (d-10)  and  (d-11) 


(d-18) 
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APPENDIX  E 

SUBROUTINE  ZIN 

SUBROUTINE  ZIN  is  used  to  compute  the  input  impedance  Z  of  a  probe 
inserted  into  a  TEM  cell.  ZIN  calls  upon  four  major  subroutines.  These 
are  (1)  SUBROUTINE  ZC  which  computes  the  characteristic  impedance  of  a 
section  of  RCTL  according  to  (18);   (2)  SUBROUTINE  HMSUM  which  sums 
hm(a,d),  and  hm(a,d)Jo[t(a2  +  M^] ;   (3)  SUBROUTINE  XSUM  which  finds  the 
reactance  X  given  by  (20),  and  (4)  SUBROUTINE  INTEGR  which  integrates 
a  complex  function  over  the  real  axis.  INTEGR  is  called  twice 
in  order  to  evaluate  the  principal  value  integral  defined  by  (21).  The 
input  to  ZIN  is  mostly  via  a  common  block  PARAM  which  contains  the  cell 
dimensions  a-,,  a?,  b.,  ,  b?,  and  g,  the  probe  dimensions  d  and  t,  the 
constant  it  ,  and  the  TEM  mode  propagation  constant  k.  Figure  6  gives  a 
flow  chart  of  SUBROUTINE  ZIN.  The  only  direct  input  to  ZIN  is  e     which 
sets  the  desired  accuracy.  ZIN  outputs  Z  ,  R,  X,  and  6X  as  well  as  the 
total  impedance  Z. 

The  constant  e  is  used  to  set  the  upper  summation  limit  in  the 
various  sums  evaluated  by  the  subroutines.  The  infinite  integral  in  6X 
is  also  truncated  according  to  e.  There  are  three  basic  sums;  the  sums 

inherent  in  L(a)  and  X,  and  the  summation  of  the  h^c^d).  We  will 

2    ?    2 
begin  with  the  last.  For  large  m,  that  is  Mn  »  kc  -   a  ,  we  have 

(1) 
KUMM, ,  and 

sin  °^  cosh  NLd    sin  °^  -Mjb.-d) 

h(a,d)^— 5 ^ l—   * ^-e  '   '  (e-1) 

m       M^  sin  h  M]b]       M^ 

We  wish  to  truncate  the  summation  so  that  the  remainder  is  less  than  e  . 
Thus  we  need  to  look  at  the  remainder  sum 

-   4a?(-l)m    -(2m-l)Tr(b1-d)/2a1 

I       — ! o-2  e  <e  (e"2) 

m=K  (2m-irTT 

Clearly  the  sum  is  dominated  by  the  following  much  simpler  sum. 
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n  /  \    v  r   t\M     -(2m-l)x 
Mx)  =  L     H)  e 
K     m=K 

so  we  will  set  M  by  requiring  that 

4a?      Tr(b,-d) 

K   ( o!—  )  <  e 


2       K  v   2a, 

IT  I 


(e-3) 


(e-4) 


Summing  R,.(x)  is  straightforward  and  yields 


RK(x) 


2e-2(K-l)x 


cosh  x 

Thus  inserting  this  result  into  (e-4)  and  solving  for  K  we  find 

2 


K  > 


1 
rr(brd) 


In 


8a 


ett  cosh 


(b~-d) 

~2lu 


+  1 


(e-5) 


(e-6) 


This  condition  will  be  used  to  truncate  the  h  (a,d)  summation. 

m 

Next  consider  the  summation  in  L(a)  given  by 

C( 


-    I 

J   mo 


:ot  K^    -,  1 
m  +  -L  f 


K 


m 


J  . 


(e-7) 


For  large  m  the  remainder  will  behave  according  to 


2   I 


*  m=K.  (2nl"1) 


1  -  coth 


(2m-l)  rrb 


\ 


2a 


(e-8) 


Again  we  choose  a  simpler  dominating  sum.  For  large  m,  such  that 
(2m-l)  rrb./2aj  >>  1,  the  coth  term  may  be  replaced  by 


coth 


(2m-lH  b 


J 


2a. 


=  1  +  2e 


-(2m-l)7rb./a. 

J    J 


(e-9) 


Therefore  we  need  to  consider 


4  \        _i_  -(Zm-lhbj/aj 


tt  m=K.  2m-l 


(e-10) 
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Clearly  we  have 


4     °°  -,        -(2m-l)TTb./a.       4       «       -(2m-l)^b  /a 

-    I       ~— -y  e  JJ<-£e  JJ 

m=K.  m=K. 


(e-11) 


Summing  the  R.H.S.  and  requiring  that  the  result  is  less  than  e/2  yields 


a  . 


K.    > 


j  27Tb. 


ix\ 


16 


-  £tt  sin  h 


rb. 


J 


+    1 


(e-12) 


This  condition  will  be  checked  for  both  j  =1  and  2,  and  the  maximum  value 

K.  chosen. 
J 

We  next  consider  the  series  contained  in  X.  Because  both  kd  and  kt 

2 
are  small,  the  decay  of  the  series  will  be  primarily  due  to  the  N-,  in  the 

denominator.  Thus,  we  will  consider  the  remainder  series  dominated  by 


I 


1 


n=N  N 


(N1  =  niT/b.,) 


(e-13) 


This  series  in  turn  may  be  bounded  by  a  simple  integral 

n=N   N  J   x 

n  im  i^        N_-,  x 

So  requiring  that  the  R.H.S.  be  bounded  by  e  implies  that 
L2 


1 


£TT 


+  1 


(e-14) 


(e-15) 


What  remains  is  to  bound  the  truncation  of  the  principal  value  integral 

-2 
In  terms  of  a   we  see  h  (a,d)  behaving  roughtly  as  a  .  Therefore,  the 

111  /- 

total  integrand  behaves  approximately  as  a~  .  Thus,  the  truncation  gives  a 
remainder  integral  similar  to 


-6  ,    C  -5 

a   da  =  -p  x 
o 

2   2     2 
where  C  =  -couk  esc  kd/4a-,  .   If  we  bound  this  result  by  e/2,  x  must 

satisfy 


(e-16) 


x  >  (2C/5e) 


1/5 


(e-17) 
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which  is  the  criterion  used  to  set  the  upper  limit  of  integration. 

A  listing  of  the  subroutines  follows.  They  were  run  at  the  University 
of  Colorado  using  the  MNF  compiler,  Version  5.4.  It  should  be  stressed  that 
e  sets  a  percentage  error  rather  than  an  absolute  error.  That  is  £  =  0.01 
implies  1%  error  rather  than  accuracy  to  two  digits  beyond  the  decimal  point, 
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Figure  3. 


Radiation  resistance  for  a]  =  a2  =  b,  =  b~  =  1  m, 
t  =  1  mm,  and  f  =  1  MHz 
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Figure  4.  Reactance  for  a.,  =  a2  =  b.  =  b~  =  1  m,  t  =  1  mm,  and  f  =  1  MHz 
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Figure  5.  Perturbation  reactance  for  a^  =  a2  =  b-j  =  b^  =  1  m, 
t  =  1  mm,  and  f  =  1  MHz 
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final  reports  on  work  performed  by  NBS  for  outside  sponsors 
(both  government  and  non-government).  In  general,  initial  dis- 
tribution is  handled  by  the  sponsor;  public  distribution  is  by  the 
National  Technical  Information  Services,  Springfield,  VA  22161, 
in  paper  copy  or  microfiche  form. 
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